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We evaluate the strangeness-conserving NN, EE, EH, AE and the strangeness-changing AN, EJV, 
AH, EH axial charges in lattice QCD with two flavors of dynamical quarks and extend our previous 
work on pseudoscalar-meson-octet-baryon coupling constants so as to include 7rHH, KAE, and i\EH 
coupling constants. We find that the axial charges have rather weak quark mass dependence and 
the breaking in SU(3)-flavor symmetry is small at each quark-mass point we consider. 
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I. INTRODUCTION 

Hyperon axial charges are significant parameters for 
low-energy effective description of baryon sector as they 
enter in the loop graphs of chiral perturbation the- 
ory. While the nucleon axial charge can be precisely 
determined from nuclear /3-decay (the modern value is 
9A,NN = 1.2694(28) 0), we do not have enough informa- 
tion about hyperon axial charges from experiment. The 
theoretical estimates from chiral perturbation theory 0- 
01, large N c limit [| of QCD and QCD sum rules [|[ exist 
but the results from these approaches are rather impre- 
cise. The lattice calculations of the axial charge of the 
nucleon have reached a mature level 043 however only 
recently there have been attempts to extract the hyperon 
axial charges using lattice QCD [Tol. ITT1]. 

In the SU(3)-flavor [SU(3)f] symmetric limit, one can 
classify the axial charges of baryons in terms of the con- 
stants of two types of couplings, F and D [l2|, as follows: 

9a,nn = F + D, 9a,se =D-F, g A ,ss = 2F, 
9A,m = 3F - D, g A ^~=-{F + D), 

9A,AN = 3F + D, 9A,T,N = D - F, 9A,AT, = 2D. 

(1) 

This systematic classification, which phenomcnologically 
works rather well but is not known a priori to hold, is 
expected to govern all the axial charges however as we 
move from the symmetric case to the realistic one, the 
SU(3)f breaking occurs as a result of the s-quark mass. 
The broken symmetry no longer provides a pattern for 
the couplings, and therefore they should be individually 
calculated based on the underlying theory, QCD. 

Recently we have extracted the xNN, 7rEE, 7rAE, 
KAN and K'SN coupling constants and the correspond- 
ing monopole masses in lattice QCD with two flavors 
of dynamical quarks [l3|. We have found that the 
SU(3) f parameters have weak quark-mass dependence 
and thus the SU(3) f symmetry is broken by only a few 
percent. Our aim in this work is two-fold: We first 



concentrate on the coupling constants 7rSS, if AS and 
iiES in order to complete our program of calculating 
the pscudoscalar-meson-octet-baryon coupling constants 
from lattice QCD that we started in Ref. [13|. In the 
second part we evaluate the strangeness-conserving NN, 
EE, SS, AE and the strangeness-changing AJV, EJV, AS, 
ES axial charges in lattice QCD with two flavors of dy- 
namical quarks. The evaluation of the coupling constants 
and the axial charges allows us to check whether the 
SU(3)f relations arc well respected in the degenerate 
quark-mass case and to what extent this symmetry is 
broken as we restore the physical masses of quarks. We 
assume exact flavor-SU(2) symmetry and take u and the 
d quarks degenerate. 



II. THE FORMULATION AND THE LATTICE 
SIMULATIONS 



We refer the reader to Ref. [T3[ for the lattice formula- 
tion and the details of the calculations of pseudoscalar- 
meson-octet-baryon coupling constants. As for the axial 
charges, we consider the baryon matrix elements of the 



isovector axial-vector current 
which can be written in the form 



(B(p)\A^\B'( P ')) =C B& u{p) [w 5 G A>BBI (q 2 ) 



+75" 



m B + mj3> 



■Gp,bb' (q 2 



u(p) 



(2) 



where 



Pu 



p^ is the transferred four-momentum 



and u(p) denotes the Dirac spinor for the baryon with 
four-momentum p and mg. Ga,bb'(p 2 ) an d Gp^B'iP 2 ) 
are the baryon axial and induced pseudoscalar form fac- 
tors, respectively. The isospin factors Cbb' are given 
as C NN = C ss = C7v S = C ss = 1, C ss = — 1/V2, 
Cae = Cah = —Cna = 1/ \/6. 

The baryon axial charges are defined as the axial 
form factors at zero-momentum transfer, viz. 9a,bb' = 
Ga,bb'{^). We compute the matrix element in Eq. @ 
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using the ratio 
R(t2,ti;p',p;T-fj,) = 

(J*A,B'(t 2 ,t i; p',p;r)) 

(FB'(i 2 ;p';r 4 )) 



(f B (t 2 -t i;P ;r 4 )) 

(^'(t 2 -t i;P ';r 4 )) (3) 



(F B '(t i; p';r 4 ))(F B '(t 2 - 7 p';r 4 )) 
(FB(t i; p;r 4 ))(Fe(t 2 ;p;r 4 )) 



1/2 



where the baryonic two- and three-point correlation func- 
tions arc respectively defined as 



(F 6 (t;p ; r 4 ))=^ e ^ x rr' 

X 

x (vac|Tfog(z)7$(0)]|TOC>, 
(F BA » B '(t 2 ,t 1 ;p',p;T)) = -i £ e ^ 

x 2 ,xi 

xr QQ '(vac|T[^(x 2 )A A1 (c Cl )^:(0)]|vac), 



(4) 



(5) 



with r = 7375T4 and r 4 = (1 + 7 4 )/2. The baryon 
interpolating fields are given as 



Vn(x) 


_ ^a-bc 


[u Ta ( 


mix) 


_ ^a-bc 


[s Ta { 




_ ^a-bc 


[s Ta { 


Va(x) 


l 

= Te 





u ra (x)C-f 5 d b (x)]u c (x), 
(x)C l5 d b (x)] S c (x), 
(x)C l5 u b (x)]u c (x), 



(G) 



{{u la {x)Cj 5 s b {x)]d c {x) - [d la {x)C 



x l5 s b {x)]u c (x) + 2[u Ta (x)C-f 5 d b (x)]s c (x)}, 

where C = 7 4 7 2 and a, b, c are the color indices. t\ is 
the time when the meson interacts with a quark and i 2 is 
the time when the final baryon state is annihilated. The 
ratio in Eq. ([3]) reduces to the desired form when t 2 — t\ 
and ij » a, viz. 



R(t 2 ,h;0,p;T;n) 



ti»q 
*2— ti^a 



E 



2m 



G 



A,BB 



(Q 2 ), (7) 



where m and E are the mass and the energy of the initial 



baryon and Q 



-q . We apply a procedure of seeking 



plateau regions as a function of ii in the ratio (O and 
calculating the axial form factors Ga,bb'(Q 2 ) at Q 2 = 
in order to extract the axial charges Qa,bB' ■ 

We employ the same lattice configuration as in our 
previous work in Ref [Hj]. It is a 16 3 x 32 lattice with 
two flavors of dynamical quarks and the gauge configura- 
tions we use have been generated by the CP-PACS col- 
laboration (l4j with the renormalization group improved 
gauge action and the mean-field improved clover quark 
action. We use the gauge configurations at f3 = 1.95 with 
the clover coefficient csw = 1-530, which give a lattice 
spacing of a = 0.1555(17) fm (a" 1 = 1.267 GcV) as de- 
termined from the p-meson mass. The simulations are 
carried out with four different hopping parameters for 
the sea and the u,d valence quarks, ft sea , «;"' , = 0.1375, 
0.1390, 0.1400 and 0.1410, which correspond to quark 



masses of ~ 150, 100, 65, and 35 McV, and we use 490, 
680, 680 and 490 such gauge configurations, respectively. 
The hopping parameter for the s valence quark is fixed to 
K vai = 0.1393 so that the Kaon mass is reproduced [l4[, 
which corresponds to a quark mass of ~ 90 MeV. We 
employ smeared source and smeared sink, which are sep- 
arated by 8 lattice units in the temporal direction. Source 
and sink operators are smeared in a gauge- invariant man- 
ner with the root mean square radius of 0.6 fm. All the 
statistical errors are estimated via the jackknife analysis. 
The renormalization factors relevant to the axial currents 
are all computed in a perturbative manner: Za =0.2576, 
0.2491, 0.2434 and 0.2377 at re = 0.1375, 0.1390, 0.1400 
and 0.1410, respectively [lij . 

Here we mention the systematic errors that could enter 
our results. Possible systematic errors arise from (1) the 
finite volume, (2) the perturbative estimation of renor- 
malization factors, (3) the finite lattice spacing, (4) the 
wrong number of dynamical quarks, and (5) the unreal- 
istic heavy quarks. 

For the finite volume effect, the present spatial lattice 
extent is 16 (about 2.5 fm) and the pion mass ranges 
from 0.440 to 0.899 in lattice unit (from 550 MeV to 
1.15 GeV), which gives 7 < m^L < 14. Serious finite 
volume artifact for the nucleon axial coupling seems to 
appear only when m^L < 7 (l5| : so we expect small finite 
volume effects in our present calculations. 

We estimate the renormalization factors in a pertur- 
bative way, which gives rise to O(10)% errors in the 
nucleon vector charge or the pseudoscalar meson de- 
cay constants [l6|. In order to reduce such systematic 
errors, we evaluate the ratios of the couplings, which 
would be less dependent on these factors. For consis- 
tency check, we compare our results with the those of 
Lin and Orginos [l(J and find very good agreement (See 
the discussion below). Then, the systematic errors from 
(l)-(4) are considered to be small in comparison with 
statistical errors. 

Although our (light) quarks are heavier than real u,d- 
quarks, we are mainly focused on SU(3) f breaking pat- 
tern. Our present setup covers a wide range of quark 
masses including the SU(3) f symmetric point (m u = 
uid = m s ). On the other hand, it is difficult to estimate 
the breaking pattern at the chiral limit. The physics at 
the physical point could be accessed with much lighter 
quarks as well as more realistic lattice setups, which is 
planned for a future work. 



III. RESULTS AND DISCUSSION 

We first concentrate on the coupling constants g^ss, 
9k as and gxT.'s- The SU(3)f relations for these pseu- 
doscalar coupling constants are given as 



-g(l - 2a), g K M 



V3 



g{4a - 1), g K ss 



-9, 
(8) 
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TABLE I: The fitted values of the 7rH3, KAz, and _KEH coupling constants and the corresponding monopole masses normalized 



with g wNN 


and A^jviv, respectively. Here, we 


define gf IBB , = g M BB> 1 


5ttJvjv and A MBB , 










9k AS 


9ksb 


A R — 


A KAS 




0.1375 
0.1390 
0.1393 
0.1400 
0.1410 


-0.227(18) 
-0.216(14) 
-0.217(14) 
-0.245(13) 
-0.273(26) 


0.334(15) 
0.348(16) 
0.347(16) 
0.313(14) 
0.291(25) 


-1.025(20) 
-1.037(18) 
-1.036(19) 
-0.998(10) 
-0.963(48) 


0.687(130) 
0.889(206) 
0.882(208) 
0.825(111) 
0.896(148) 


1.030(164) 
0.908(149) 
0.918(157) 
1.085(146) 
1.186(242) 


0.884(39) 
0.891(44) 
0.891(49) 
1.044(28) 
1.237(99) 
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FIG. 1: The NN axial charge together with S3, EE, AE 
axial charges normalized with qa,nn as a function of mj. 
The empty circle denotes the SU(3)f limit. 



FIG. 2: Same as Fig.[T]but for strangeness-changing AH, EH, 
AN and EA axial charges. 



where g = g^N n and a is the F j (F + D) ratio of the 
pscudoscalar octet. We extract these coupling constants, 
Smbbs an d the corresponding monopole masses, A MBB >, 
for each ft™';. Our results are presented in Table HJ We 
give the fitted values of the 7rSS, KAE, and KT,E cou- 



pling constants and the corresponding monopole masses 
normalized with g^NN and A w jvjv ; respectively. In 
Table [fl gf IBB , and Af IBB , denote gMmi/dnNN and 
A A/eg' /Attjvjv, respectively. As in Ref. [13| we expect 
that the systematic errors cancel out to some degree in 
the ratios of the coupling constants and those of the 
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monopolc masses. 

In the SU(3)_f limit, where K^' a[ = n s val = 0.1393, the 
SU(3)f relations in Eq.® are exact and all the cou- 
pling constants are well reproduced with a = 0.395(6), 
which is obtained by a global fit including the coupling 
constants obtained in Rcf. In the SU(3) f broken 

case, we observe that our conclusion in Ref. [13| for the 
pseudoscalar-meson-baryon coupling constants holds as 
well for the coupling constants in question here: The 
quark-mass dependences for gjy-gg/ and Aj^g B , are not 
large and the ratios of the coupling constants, <?mbs'' are 
similar in value to those in the SU(3)f symmetric limit, 
and the monopole-mass ratios, A^ /BB , , are almost unity 
independently of the quark masses. This confirms that 
the SU(3)f breaking is small for pseudoscalar-meson- 
baryon coupling constants at the quark masses we con- 
sider. 

We next concentrate on the strangeness-conserving 
and strangeness-changing axial charges of baryons. The 
axial charge of nucleon in the present setup is about 10% 
overestimated Q , which would be due to the perturbative 
estimation of renormalization factors 

[Hj]. The present 

lattice spacing is about 0.15 fm, which is far from the con- 
tinuum limit. In fact, the vector charge of nucleon as well 
as the decay constants obtained with the same setup as 
ours show O(10)% deviation from physical values fl4lfl6l ]. 
Hence, we evaluate the ratios of axial charges (charges 
normalized with the axial charge of nucleon) rather than 
the bare values, so that we expect the cut-off artifacts in 
the renormalization factors to cancel to some extent. 



K vai = <ai = 0.1393, we obtain a A = F/F + D = 
0.390(2). The value of a a has a weak quark-mass depen- 
dence and as we approach the chiral point a a tends to 
decrease. We illustrate this behavior in Fig. [3] 



.1 1 1 1 1 1 1 1 1 1 1 1 1 




1 1 1 


1 | 1 1 1 | 1 1 




r I 3 5 a 
- j 1* 




* * 


* 









B = 


E [Linetal.] 






• 


B = 


3 [This work] 






A 


B = 


£ [Linetal.] 






A 


B = 


£ [This work] 




- g D *a □ » 




* • 






:, , , 1 , , , 1 , , , 1 , 




1 1 1 


■ ■ . . . 1 . . 





0.2 0.4 0.6 0.8 1 1.2 1.4 

ml [GeV 2 ] 



FIG. 4: Com par ison of our results for g^,EE an d 5a,hh with 
those in Ref. [10f | 

In Fig. HI we compare our results for g R ss and g R ss 
with those obtained from a lattice setup using stag- 
gered fermion action for the sea quarks and domain- wall 
fermions for the valence quarks [10]. We observe the re- 
sults from two different setups are in good agreement 
with each other. 
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FIG. 3: a.A = F/F + D ratio as a function of m\. The empty 
circle denotes the SU(3)f limit. 

In Table HH we give the fitted values of the NN ax- 
ial charge, g A ,NNi together with the fitted values of 
the strangeness-conserving SS, EE, AE and strangeness- 
changing AS, ES, AN and EiV axial charges normalized 
with gA.NN for various quark masses and illustrate these 
in Figs. [1] and Similarly to the pseudoscalar-meson- 
baryon coupling constants in our previous analysis [13 ]. 
we expect that the systematic errors cancel out to some 
degree in the ratios of the axial charges. Here, we define 
9a bb 1 ~ 9a,bb' I 9a,nn ■ We also present the values of 
the ratios of the coupling constants a a = F/F + D as 
obtained from a global fit. In the SU(3)f limit, where 



FIG. 5: The SU(3)-breaking parameter, <5su(3), as a function 
of m^. The empty circle denotes the SU(3)_f limit. 

In the SU(3)f broken case, the deviations in the cou- 
pling constants are not large and the values are similar 
in those at k = 0.1393. In order to quantify the SU(3)i? 
breaking, we construct the following combinations: 





- 9a,ss + 


9a,t,si A 2 -- 


= 2$£sB + flA,AE J 


-h 


= (9%ss 4 


- 5a,ajv)/2, 


Ai = 9a, as - 5a,hh> 


A 5 


= ^9a,ts ' 


~ 9a,AE\i A 6 


- n R n R 

— 9 A, AN #A,££) 


A? 


= 9a,sn + 




= (9a,as+3a,se)/2, 


A 9 


= (^9a,an 




A W = 9A.AE, + 2 9a,XN> 



5 



TABLE II: The fitted value of the NN axial charge together with the fitted values of the strangeness-conserving S3, EE, AE 



and strangeness-changing AH, EH, AN and EiV axial charges normalized with Qa,nn • 
We also give the fitted value of F/F + D at each quark mass. 



Here, we define g ABB , = g A ,BB> 



9a,nn 



ah 



n R 



n R 

Pa. as 



3.4.AH 



3a.se 



9a,an 



n R 

9 A, SJV 



F/F + D 



0.1375 
0.1390 
0.1393 
0.1400 
0.1410 



1.284(11) 
1.282(15) 
1.280(15) 
1.289(15) 
1.314(24) 



0.218(05) 
0.220(04) 
0.221(04) 
0.221(04) 
0.228(06) 



0.791(04) 
0.782(04) 
0.779(04) 
0.772(04) 
0.738(09) 



1.223(05) 
1.221(04) 
1.221(04) 
1.218(04) 
1.221(12) 



0.564(09) 
0.558(08) 
0.559(09) 
0.553(07) 
0.511(14) 



-0.994(03) 
-0.999(01) 
-1.000(01) 
-1.000(02) 
-0.977(11) 



1.761(06) 
1.776(04) 
1.779(04) 
1.790(05) 
1.775(14) 



0.212(04) 
0.220(04) 
0.221(04) 
0.225(04) 
0.258(08) 



0.390(2) 
0.390(2) 
0.390(2) 
0.390(2) 
0.380(3) 



An = (sIah + 2<?I AS )/3, A 12 = (g* + g^ N )/2, 



^13 - (9a,an 



Air 



9a,at,)/3, A u = g A A s - g% 



In the SU(3)f symmetric limit, the above equations sat- 
isfy Ai = A<2, = . . . ee A15 = 1, which can be verified 



by inserting the coupling constants at k 1 z 



0.1393 in 



Table UU At other quark masses, the deviations from 
unity represent the amount of SU(3)f breaking. Insert- 
ing the values of the coupling constants corresponding to 
the lowest quark mass we consider in Table U into (jH]) , we 
find Ai = 0.966(07), A 2 = 0.966(11), A 3 = 1.002(07), 
A 4 = 0.993(08), A 5 = 0.965(10), A 6 = 1.037(10), 
A 7 = 0.996(05), A 8 = 0.979(06), A 9 = 1.013(07), 
Aio = 1.027(10), A n = 0.984(08), A 12 = 1.017(06), 
A13 = 0.999(06), A 14 = 0.962(10), A 15 = 0.977(10) 
which indicate a breaking in SU(3)_f by less than 10%, 
as we approach the chiral limit. Moreover, we define the 
average SU(3) f breaking as follows: 



3SU(3) 



A n \ 



(10) 



which amounts to <5 SU(3) =0.010(1), 0.002(1), 0.006(1), 
and 0.021(4) for the quark masses at ~ 150, 100, 65, and 
35 MeV, respectively. This suggests for the axial charges 
of the octet baryons that SU(3)f is a good symmetry 
in the quark-mass range we consider, which is broken by 
only a few percent, similarly to the pseudoscalar- meson 
coupling constants. We have also tried a quadratic fit 
°f ^su(3) an d extracted #su(3) = 0.047(3) in the chi- 
ral limit. Fig. [5] shows the value of <5su(3) as a func- 
tion of m^. and the chiral extrapolation with errors. As 
for a a = F/F + D, it seems to have a slightly smaller 
quark-mass dependence as compared to a = F/F + D of 
pseudoscalar-meson-baryon coupling constants [l3j . It 
is interesting to note that cha as extracted from axial 
charges is closer to 2/5, prediction from SU(6) quark 
model, in the present quark-mass range. 

Assuming pion-pole dominance, the axial charges of 
octet baryons are related to their pseudoscalar-meson 
coupling constants via Goldberger-Treiman relations: 



fM9MBB' = SBB'{m B + ■m B ,)/2g A ,BB> 



a 1.1 r 




Tfln |lattice unit| 



FIG. 6: Jbb' as defined in Eq. (fT3")l as a function of : 



where is the meson decay constant and Sbb' ar e fac- 
tors that fix our convention in choosing the isospin fac- 
tors in Eq.Q with respect to Goldberger-Treiman re- 
lations: SnN = EE 6>£AT EE 5eh EE — 6>EE = 1, 
Sad ee — San = Sae = l/v3. We can check the va- 
lidity of these relations by normalizing both sides with 
the Goldberger-Treiman relation of the nucleon and pion, 



firg-rrNN = mjy gA,NN, 

and by constructing the following form: 



Jbb 1 — Si 



BB' 



9a,bb> U (rag + m B ,) 
9mbb' f M 2mN 



(12) 



(13) 



in order to eliminate the systematic errors. If these rela- 
tions are valid, then we should have Jbb 1 =1- In Fig. |6l 
we present Jbb 1 for all the couplings in question as a 



function of 



Since we have large error bars, it is 



(11) 



difficult to reach a definite conclusion about the valid- 
ity of the Goldberger-Treiman relations, however, we ob- 
serve that these relations are rather good in the SU(3)p- 
symmetric limit and discrepancies arise, in particular for 
5 couplings, as wc approach the chiral limit. 
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IV. CONCLUSION 

We have evaluated the strangeness-conserving NN, 
EE, SS, AE and the strangeness-changing AN, UN, AS, 
ES axial charges in lattice QCD with two flavors of dy- 
namical quarks and extended our earlier work in Ref. [l3f 
so as to include the 7rSS, if AS and ifES coupling con- 
stants. We have extracted the ratios of the axial charges, 
which are supposed to be less prone to systematic errors. 
The ratios of the axial charges show rather weak quark- 
mass dependence. We have allowed an SU(3)_f breaking 
by varying the quark masses. Our results suggest that 
the SU(3) f for axial charges is a good symmetry, which 
is broken by only a few percent. This conclusion is in 
agreement with what we have found for the pseudoscalar- 
meson couplings of the octet baryons. While we think 
that the present work reveals the SU(3) f pattern of axial- 
vector couplings of octet baryons, there are a number of 
improvements to be considered in a future work. Our lat- 
tice is still coarse by modern standards and quark masses 
are too large to reach a definite conclusion about SU(3) f 



breaking. Simulations with more realistic setups with 
smaller lattice spacing and larger lattice size employing 
muc h lig hter quarks and a dynamical s-quark arc under 
way [17] . It is an intriguing issue to see whether our find- 
ings in this work are retained in more realistic 2+1-fl.avor 
calculations with much lighter quark masses. 
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